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Abstract 

The theoretical basis of the phenomenon of effective and exact dimensional 
reduction, or holographic correspondence, is investigated in a wide variety 
of physical systems. We first derive general inequalities linking quantum 
systems of different spatial (or spatio-temporal) dimensionality, thus estab- 
lishing bounds on arbitrary correlation functions. These bounds enforce an 
effective dimensional reduction and become most potent in the presence of 
certain symmetries. Exact dimensional reduction can stem from a duality 
that (i) follows from properties of the local density of states, and/or (ii) from 
properties of Hamiltonian-dependent algebras of interactions. Dualities of 
the first type (i) are illustrated with large-n vector theories whose local den- 
sity of states may remain invariant under transformations that change the 
dimension. We argue that a broad class of examples of dimensional reduction 
may be understood in terms of the functional dependence of observables on 
the local density of states. Dualities of the second type (ii) are obtained via 
bond algebras, a recently developed algebraic tool. We apply this technique 
to systems displaying topological quantum order, and also discuss the impli- 
cations of dimensional reduction for the storage of quantum information. 
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1. Introduction 

Spatial (or space-time) dimensionality, topology, and symmetry are cen- 
tral concepts in physics. This article studies a phenomenon, that of dimen- 
sional reduction or holographic correspondences, that highlights how these 
concepts may come together to produce particularly striking signatures. 

The dimension of a system is one of its most basic characteristics and 
typically represents an important measure of its complexity. Dimensional 
reduction and holographic correspondences refer to the fact that the "ap- 
parent", or "obvious" dimension D of a system may not be the dimension 
d, d < D, that best characterizes its response to experimental probes and 
its information content. We speak of dimensional reduction when the ob- 
servable physical properties of a system in D dimensions behave as if it were 
d-dimensional. We argue for a holographic correspondence when two systems 
of different dimensionality seem capable of storing and processing compara- 
ble amounts of information. Both notions, often undistinguished in current 
literature, appear in numerous fields including condensed matter, cold atom, 
high-energy, and black-hole physics. Given the cross-disciplinary nature of 
the subject there is a need to better understand the physical and mathe- 
matical basis of dimensional reduction and holographic correspondences. We 
believe that this understanding should be rooted in basic principles than in 
properties of specific models. 

This article constitutes an attempt to study exact and effective dimen- 
sional reductions, and holographic correspondences within a unified frame- 
work. Our approach relies on the recent theory of dualities [1], and a new 
theorem that constrains arbitrary correlators of a quantum theory in D di- 
mensions. The obtained inequalities and bounds afford a practical notion 
of effective dimensional reduction because they are specified by homologous 
correlators of an associated, local, d < D dimensional theory. Also they be- 
come specially useful in the presence of d- dimensional gauge-like symmetries 
(or d-GLSs for short) [2, 3], with d < D, a geometrically distinguished set of 
symmetries that we will discuss further below. The notion of exact dimen- 
sional reduction is further linked to the existence of a duality, mapping the 
-D- dimensional theory to a d- dimensional one. In such situations inequalities 
may be replaced by exact equalities. We also introduce the concept of di- 
mensional reduction by the density of states, a concept that may find its way 
to practical applications in the form of either exact or effective versions, de- 
pending on the problem. We illustrate this technique with 0(n) type vector 
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theories. Finally, we will show that exact dimensional reduction is sometimes 
connected to a holographic-like entropy that scales as the surface of the sys- 
tem, and use large- n vector theories to illustrate this point. Remarkably, our 
take on dimensional reduction/holographic correspondences will allow us to 
discuss limitations on the storage and robustness of quantum information. 

While dimensional reduction/holographic correspondences may seem to 
be exotic phenomena, they can in fact appear in many systems as the natural 
consequence of a wealth of physical mechanisms: 

(1) Restrictions from conservation laws. In some cases, a conservation 
law (that can be formalized as a gauge constraint) favors sliding dynamics 
along preferred directions while forbidding motion along other directions. 
For example, the glide principle [4] for elastic solids dictates that dislocations 
move preferentially along directions that are determined by their topological 
charges or Burger's vectors. 

(2) Reduced kinetics and interactions along one or more directions. In 
many condensed matter systems a dimensional reduction can occur due to the 
confinement of electrons/atoms along planes or lines (e.g., wires), whereby 
certain directions become irrelevant. Reduced tunneling, for instance, can 
appear along certain directions, as in the case of high-T c cuprates and some 
other layered systems, due to relatively large separation between planes that 
suppresses interplane couplings and tunneling from one layer to the next. 
Similar decoupling can occur in magnetic systems with spin exchange inter- 
actions on geometrically frustrated lattices [5], where the coupling between 
planes (or other lower- dimensional subvolumes) along one or more directions 
becomes negligible. In all of these systems, the effective physical interactions 
and/or kinetic hopping matrix elements describe a lower- dimensional system. 
Many of these systems are locally inhomogeneous with electrons constrained 
in some way, but at large distances the symmetry between different directions 
is restored. That is, in such cases, converse to the Kaluza-Klein type com- 
pactifications that we list next, dimensional reduction may (but need not) 
emerge at local scales. At short distance scales the system is effectively of a 
lower dimensionality than it is at large distance scales. 

(3) Kaluza-Klein compactification. Dimensional reduction can be achieved 
geometrically by compactifying (D—d) spatial dimensions of a theory through 
the Kaluza-Klein [6] procedure. At low energies or temperatures, excitations 
along these "small dimensions" are suppressed, leading to an effective d- 
dimensional theory. That is, dimensional reduction emerges at low energies 
and large distances. This type of dimensional reduction is important in string 
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[7] and supergravity theories, and it has been suggested recently [8, 9] that 
may also help to address the Langlands duality which is an open problem in 
basic mathematics. 

(4) Gravity-gauge dualities (AdS-CFT). A seemingly different, more sub- 
tle, and very powerful form of dimensional reduction is offered by the AdS- 
CFT (or more general gravity-gauge) holographic correspondences (duali- 
ties). This correspondence relates, via a putative strong- weak coupling du- 
ality, different theories: string theory as a theory of quantum gravity on an 
Anti de Sitter (AdS) background in D dimensions, and a conformal field 
theory (CFT) in d = D — 1 space-time dimensions. Since its initial discov- 
ery [10], the AdS-CFT correspondence [10, 11, 12, 13, 14] has been a source 
of inspiration for many ideas and applications. Apart from its original use 
in string theory, as a gravity-gauge correspondence [10], it has led to many 
other applications that include hadronic hydrodynamics [15] and fundamen- 
tal questions relating to transport in strongly correlated systems [16]. Similar 
forms of an AdS-CFT type correspondence were conjectured for cold atom 
systems with Schrodinger type dynamics [17]. 

This incomplete list of known mechanisms underscores the multitude of 
disparate phenomena that may lead to dimensional reduction and holographic- 
like correspondences. This article will describe simple tools that can help ad- 
dress this diversity in a coherent way. As mentioned already, we distinguish 
between two scenarios: 

(i) Effective dimensional reduction. This possibility has, for the most 
part, not been explored before (especially in quantum theories), and forms a 
cornerstone of this article. Items (1), (2), and (3) above can be seen as real- 
izations of this type of dimensional reduction. We will illustrate how inequal- 
ities with bounds on arbitrary correlation functions, typically strengthened 
by the presence of d-GLSs, relate to an effective dimensional reduction. 

(ii) Exact dimensional reduction. We identify this scenario with the ex- 
istence of dualities connecting theories of different dimensionality. These 
dualities may be exact or asymptotically so in some limit (such as the large- 
n limit or the weak coupling limit). Case (4) is often conjectured to be 
exact yet still largely remains to be proved. The AdS-CFT correspondence 
has, however, been convincingly established in the large-n limit and several 
related SU(n) matrix and other models. 

In general, an effective dimensional reduction is very different in na- 
ture from an exact dimensional reduction, whence the spectrum of a D- 
dimensional system attains a form akin to that of a theory in d < D di- 
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mensions with local interactions. In this paper, this effective scenario is 
characterized in terms of a flexible formalism that can describe a range of 
situations that go from "very weak" to "very strong" dimensional reduc- 
tion. It is essential then to recognize ingredients that will favor the strongest 
forms of effective dimensional reduction. Symmetries, and more specifically, 
the d-GLSs that we describe in the next section, are one such ingredient. 
In such cases, those symmetries mandate how the bounding d-dimensional 
correlation functions decay with distance. 

The concept of exact dimensional reduction (ii) is more restrictive, and 
its connection to d-GLS is less clearly understood. If an exact equivalence 
is present between all possible correlation functions in two different theo- 
ries then it will be natural to conjecture that a duality (i.e., an exact unitary 
mapping [1]) connects the theory that displays lower- (d) dimensional behav- 
ior, yet described as high- (D) dimensional theory, to a dual theory explicitly 
<i-dimensional. Then the inequalities saturate and become equalities. This 
occurs, for instance, in some special models [2, 18], such as Kitaev's model 
[19], that may be constructed and solved using these unitary mappings. As 
we will show, this occurs far more generally in the large-n limit of rather 
general 0(n) type vector theories. The crux of this result is that in sys- 
tems such as large-n vector theories, the partition function (or generating 
functional) of the system can be expressed in terms of an effective density 
of states functional. Transformations between systems that reside in differ- 
ent (space/space-time) dimensions that preserve this density of states, while 
maintaining the locality of the theory, automatically correspond to (dimen- 
sion reducing) duality transformations. 

We expect, though, that the most complex cases of dimensional reduc- 
tion/holographic correspondences may be understood as the combined use 
of (i) and (ii). 

2. Systems of study 

The results reported in this work hold for a very broad class of systems. 
Throughout, we will consider arbitrary quantum systems (or their classical 
descendants) in D- spatial (or space-time) dimensions with Hilbert spaces 
which are either finite or denumerably infinite. The theories to be directly 
analyzed are defined on a D-dimensional lattice of volume A. When alluded 
to, continuum field theories will be understood to be represented by the 
continuum (vanishing lattice size a) limit of an appropriately defined lattice 
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field theory regularizing it. The Hamiltonians (or Euclidean actions) of the 
systems that we address have interactions local in space (or space-time), i.e., 
have couplings that veer to zero at large spatial (or temporal) separations. In 
the thermodynamic (A — > oo) limit we will, at times, replace pertinent sums 
over the Fourier modes k associated with lattice systems by corresponding 
integrations over the first Brillouin zone. For instance, in the case of a 
square lattice system, k^ = 2ttI^/L with L the linear lattice size along the 

direction, /x = 1,2, —L/2a < Z M < L/2a being integers, with a the lattice 
constant. The first Brillouin zone is spanned by —ixja < < tt /a. 

Symmetries represent an important element of our analysis. The most 
general symmetry 14 of a system with a Hamiltonian H can be described as a 
direct sum of arbitrary unitary transformations U{NE i ) that act non-trivially 
only on their corresponding subspaces of fixed energy Ei and of dimension 

■ That is, the direct sum of unitary transformations 



constitutes a symmetry of the system, i.e., WHU = H, for arbitrary U{NE i )- 
Conversely, if U is a symmetry and H = EiPi is the spectral decomposi- 
tion of H, then U has the structure of Eq. (1) above, with U(NE i ) = UP%- 
Typically, unitary operations of the type of Eq. (1) may correspond to non- 
local operations. In contrast, the symmetries that are most relevant to di- 
mensional reduction have well defined spatial support. 

This observation (that will be illustrated often in this paper) leads to 
the notion of d-GLSs. As alluded to earlier, these are symmetry operations 
(of codimension D — d) that act non-trivially only within a d-dimensional 
physical space. Technically, a symmetry is termed d-GLS if fields acted on 
by the symmetry operators IA& are affected within a region of d-dimensional 
support [2, 3]. For instance, local (gauge) symmetries are of dimension d = 0. 
Similarly, global symmetry operators act on all of the fields of the system 
and thus their dimension is d = D (that of the entire system). In between 
these two extremes of local (d = 0) and global (d = D) symmetries there is 
an intermediate regime where the symmetry operators act on a < d < D 
dimensional spatial region. It was recently shown [2] that d-GLSs can man- 
date the appearance of topological quantum order (TQO) [20]. In continuum 
space-time d-GLSs may asymptotically include, as a particular case, the con- 
formal symmetries that appear in conformal fields theories (CFTs) but are 
not limited to these. 




(1) 
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It will be useful in the following to keep in mind specific simple examples 
for which our theorems have easy to grasp yet non-trivial consequences. The 
spin S = 1/2, D = 2 compass model [2, 3] on a square lattice is one of many 
such examples. It is often specified by the Hamiltonian 

^compass = - ^(Jl ^rc+ei + J 2 <7^ +e2 ), (2) 
aieA 

with exchange couplings J At , where <r%~ x ' y are Pauli operators located at the 
lattice site x = i x ei + i 2 e 2 (x E A). Its d = 1-GLSs 

^=n°? i ,*»' ^=ik^> ( 3 ) 

i 2 i 1 

([^compass, ^1] = = [if com pass, Yp]) satisfy the algebra 

{Xii,Y % *} = 0, [X^X?,] = [Y^Y^,] = 0. (4) 

The Xu- Moore model [21], studied before in connection to dimensional re- 
duction, is dual to the compass model [1, 22]. 

3. Effective dimensional reduction via bounds 

As explained in the introduction, our approach to effective dimensional 
reduction is based on inequalities. These inequalities are the subject of a 
theorem on an effective quantum dimensional reduction (EQDR) that we 
prove in this section. 

Intuitively speaking, a system is intrinsically D-dimensional if it cannot 
be simulated by a lower- dimensional one (a higher-dimensional simulator is 
always possible provided some redundancy is allowed [23]). This idea pro- 
vides the basis for a good definition of intrinsically D -dimensional, but fails 
to provide simple, practical criteria to search for dimensional reduction. The 
EQDR theorem provides one such criterion. The idea is that D-dimensional 
systems have observables / that specifically probe rf-dimensional sub- volumes 
r, d < D, of the Z)-dimensional space A. The EQDR theorem illustrates that 
under certain (not very stringent) conditions, the expectation values of these 
d- dimensional observables are both bounded from above and from below by 
expectation values taken within effective d- dimensional theories. Thus, we 
can use the behavior of these bounds (e.g., how stringent they are for spe- 
cific theories) as indicators to search for effective dimensional reduction. As 
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stressed in the introduction, the EQDR theorem becomes most useful in the 
presence of suitable d-GLSs. 

The EQDR theorem exploits features of quantum mechanics that force 
us to reconsider the notion of localizability. This notion poses a problem that 
was perhaps first fully appreciated in the context of entanglement [24, 25, 
26] but was also carefully considered in connection to the notion of particle 
(more specifically, the position of a particle) in quantum field theory [27]. 
Remarkably, the notions of localizability that emerged from those two lines of 
inquiry are (as far as we can see) conceptually different, and exploit different 
mathematical structures on the space of states. Hence, in order to apply 
the EQDR theorem to quantum field theories and many-body systems we 
need to introduce a redundant representation of the state space of identical 
particles, and keep the constraint of indistinguishability explicit in the form 
of a projector onto subspaces of completely (anti) symmetric state vectors. 
This enables us to treat fermions and bosons. In the current work, we do not 
discuss the extension of the EQDR theorem to systems with anyon statistics. 

3.1. A theorem on effective quantum dimensional reduction (EQDR theorem) 

The EQDR theorem is inspired by an analogous theorem for classical 
statistical mechanics [3] . Hence it is convenient to review briefly this theorem 
on classical dimensional reduction to motivate its quantum analogues and to 
introduce some useful notation. 

Consider a classical model of statistical mechanics that, for concreteness, 
has its degrees of freedom defined on the sites of a lattice A = V U A. The 
lattice provides a set of labels. In the following, we will loosely allude to V 
as "the boundary" of the system and to A as "the bulk" of the system. The 
elementary degrees of freedom of the model can be organized into a lattice 
field <p(x) (x G A denotes a site of the lattice) that can be broken up into 
two fields with disjoint supports, T fl A = 0, 



If / is a classical observable with support on V then we can write /[</>] = 
f[4>o] so that its canonical ensemble average is 




x e r 
xeA 



(5) 



e -/3E[0o,V] 



E 
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(6) 
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where E[<j>] = E[(j)Q,xjj] is the energy functional, Z = E^e - ^ 1 * 1 , with 
P = 1/ksT where ks is Boltzmann's constant, and z[ip] = X]{</> } e~^ E ^ '^. 
The quantity 

d s W» (7) 

z m 

represents a conditional expectation value of / dependent on the condition 
that within the bulk of the system A the field assumes the value ip. The 
probability distribution associated with this condition is p{ijj) = z[tp]/Z. Let 
us furthemore define 

(/)? = min^ (f) d [ij} = (/»min], (f)i = max^ (/)«ty] = (f) d [^ max \. (8) 

It then follows from Eq. (6) that [3] 

</)? < (f}° < (f)t (9) 

The key observation is that we can think of these bounds as averages for 
effective d- dimensional theories. Moreover, these theories are characterized 
by local energy functionals 

Ei[4> ,ip m \n\ and E u [<f) , Vw]- (10) 

Thus, the inequalities of Eq. (9) afford a classical notion of effective dimen- 
sional reduction. For some systems, it is, to a certain extent, a matter of 
taste whether we choose to call these inequalities "classical" or "quantum" as 
they are directly applicable to lattice quantum field theory in the Euclidean 
path-integral representation [28]. However, the path-integral formalism is 
not always convenient, and some subtleties may arise when this formalism is 
applied to quantum theories with a boundary. Thus we would like to under- 
stand next how to extend this type of reasoning to the operator formalism 
of quantum mechanics. 

Suppose that we have a quantum system occupying a D- dimensional spa- 
tial volume A. This system is described by the density matrix p that we 
take to be completely arbitrary for now (in particular, it could be time de- 
pendent). That is, the density matrix p is solely required to be a positive 
and normalized TrA(p) = 1 operator acting on the state space T-L\, with TrA 
denoting the trace on H\. A natural question to ask is whether p should 
play in the quantum arena a role similar to that played by the Boltzmann 
probability distribution in the classical setting. Another important question 
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is "what does it mean for a quantum observable / to be localized on T?". 
The answer to the latter question depends on how the set decomposition 
A = T U A of the labels is reflected in the state space 1-Lx- If, for whatever 
physical reason, the Hilbert space associated with the entire system can be 
expressed as a direct product of the spaces associated with the boundary and 
bulk, 

n A = n r ®n- A , (ii) 

i.e., if the boundary T and bulk A are distinguishable, then we can use the 
standard notion of subsystems in quantum mechanics to say that / is localized 
on T provided that 

/ = /r®lA- (12) 

Once we agree on this notion of localizability, the EQDR theorem follows 
from a simple direct computation involving partial traces. 
First we notice that 

(f}° = Tr A (p/) = Tr A Tr r (p/), (13) 

where Trp denotes the partial trace operation that takes an operator on "Ha = 
Hr ® %A to an operator on the state space for the bulk, and Tr A denotes the 
standard trace on "H A . Equation (13) is analogous to the first part of Eq. 
(6) in the classical case. To reproduce the second part, we assume that the 
reduced operator Trr(p) is invertible, so that we can write 

(f) D = ^(Trr^tTrr^)- 1 ^^/)]). (14) 

Very loosely speaking, in such a case, we can think of Trr(p) _1 Trr(p/) as a 
"conditional probability" where the conditions pertain to the quantum states 
of the bulk. Explicitly, the reduced operator Tr r (p) can be decomposed 
in terms of pure states of the bulk as Tr r (p) = £\ r ki P Ki , where P Ai is 
the projector onto the eigenspace with the eigenvalue r Ai . This (positive) 
eigenvalue represents the probability to find the bulk in the state P Ai . Then, 

(f) D = £ r * Tr A (P A ,[Tr r (p)- 1 Tr r (p/)]). (15) 

i 

To compare Eq. (15) with its counterpart Eq. (6) in the classical dimen- 
sional reduction theorem, notice that Eq. (6) expresses the classical average 
as a conditional expectation value weighted by the probabilities for different 
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possible bulk configurations. Equation (15) has a similar structure. How- 
ever, while the left-hand side of Eq. (15) is real by construction, the quan- 
tities Tr A (P Ai [Trr(p) _1 Trr(p/)]) need not be real valued since the operator 
Trr(p) _1 Trr(p/) need not, in general, be Hermitian. [This holds even if we 
assume that / is Hermitian (an assumption that we will not invoke)]. Later 
on, it will be useful to consider the situation when / is normal ([/, f>] = 0). 
Thus, we should think of Tr A (P Ai [Tr r (p)" 1 Tr r (p/)]) clS db generalized con- 
ditional average, in the spirit for instance of Dirac's work [29]. With this 
understanding we define 

(/)? = min J Re{Tr A (P A4 [Tr r (p)- 1 Tr r (p/)])} = (/)t m , (16) 
(/>; = max J Re{Tr A (P Aj [Tr r (p)" 1 Tr r (p/)])} = (/)t ax ., (17) 
Thus, from Eq. (15) it follows that 

</)? < (f}° < (f) d u ■ (18) 

Thus we have, in Eq. (18), proved the first theorem of this paper: 

Theorem 1 (EQDR). Consider a quantum system in the state p, occupying 
a P-dimensional region A, with a decompositon of A into a d- dimensional 
boundary T and its complementary bulk A, such that A = T U A. Suppose 
that the state space of the system admits a tensor product decomposition 
Ha = Ht ® H&j an d let / = fr <E> 1 A be localized on T. If the subsystem 
density matrix Trp(p) of the bulk is invertible, and if the P-dimensional 
system is local then the expectation value (f) D = Tr \(pf) can be bounded 
both from below and from above (Eq. (18)) by suitably defined expectation 
values (Eqs. (16) and (17)) of fr for local d- dimensional {d < D) systems. 

Some comments are in order. First, the range of the interactions in the 
effective theories providing the bounds is the same as in the original D- 
dimensional theory. Thus, as alluded to above, if the high (P)-dimensional 
system is local then no long-range effective interactions will appear in the 
low (<i)-dimensional system. This is different from what would generally be 
expected if a lower- dimensional description were to be achieved by brute 
force integration of bulk fields (in D dimensions) in order to generate a lower 
(d)- dimensional boundary theory. Second, we can gain some insight into 
the condition that Trp(p) be invertible by noticing, for example, that a fi- 
nite temperature density matrix corresponding to a system Hamiltonian H, 
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i.e., p = e _/3i ^/TrA(e~^ H ), is always invertible. The same applies for the 
corresponding reduced Trp(p). It is interesting to notice, however, that the 
Boltzmann distribution cannot be (easily) extended to gravitational systems 
[30] already in their classical incarnation. Third, there is a different version 
of the EQDR theorem that holds for arbitrary p, and that we expand in 
Appendix A. Fourth, while we choose to refer to A as the volume occupied 
by the system, A could, in fact, represent any suitable set of labels or modes 
without reference to any specific notion of space (or space-time). For exam- 
ple, if A denotes momentum space then our theorem will refer to localization 
and dimensional reduction in momentum space. Fifth, the EQDR theorem 
puts no restrictions on / other than / = /r <8> 1a- m particular, / does 
not have to be Hermitian. Finally, with an eye towards applications of the 
inequalities, we remark that we may choose / to be invariant with respect 
to a given group of d-GLSs {Ud}- That is, we may replace / by an operator 
/ sym which is symmetrized with respect to {Wd}, 

f^f sym = Y^U d fU d . (19) 
Wd 

3.1.1. The EQDR theorem in quantum field theory and many-body physics 

The natural state space in quantum field theory and many-body physics 
is that of Fock space which automatically takes into account the quantum 
statistics of indistinguishable particles. A Fock space can be described most 
efficiently as an irreducible representation of an algebra of operators of par- 
ticle creation (d x .) and annihilation (d x .) on a vacuum state vector. The 
actual kind of statistics is not important in what follows, so we consider op- 
erators identified by some suitable set of labels that we denote generically by 
Xi, Xj e A, and that satisfy 

ft*, <U = <W [d Xi , d x .] ± = 0, [dl, 4.] ± = 0, (20) 
where [A, B]± = AB ± BA. These operators act on a Fock space 

oo 

n Fock = U ® Hi @U 2 @--- = 0^„, (21) 

n=0 

where Ho is the subspace spanned by the vacuum state |0), and % n is spanned 
by the n-particle states 

4A---4J0} (22) 
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(aij, % = 1, • • • , n, need not be all different if we are dealing with bosons). 

The explicit description of the Fock space shows that there is a difficulty 
in applying the EQDR theorem to quantum field theory directly, since in 
this case the decomposition A = V U A of the set of labels into bulk and 
boundary is not readily reflected in the decomposition of H Fock into a tensor 
product as in Eq. (11). To find a solution to this problem, let us imagine 
for a moment that the particles in the bulk A are distinguishable from the 
particles in the boundary T, keeping however the particles in the boundary 
(bulk) indistinguishable. Then the space "Ha becomes 

H A = H F r ° ck ®H A ock = [0ftm J ® (0^An) • (23) 

\n=0 J \n=0 J 

This space has the right properties to make our EQDR theorem applicable. 
The true (irreducible) state space of the system, the Fock space K Fock , can 
be realized as a proper subspace of H\. To see this, notice that we can 
reorganize the right-hand side of Eq. (23) as 

oo 

^a = U rm ®U Am ,, (24) 

n=0 m+m'=n 

thus decomposing "Ha into subspaces with fixed numbers of particles. These 
subspaces include n-particle states that do not respect the indistiguishability 
of the particles in the whole system, but contain as a subspace 

n nC H Tm ®H Km ,, (25) 

m+m'=n 

the space of properly ( anti- ) symmetrized n-particle states. 

We can then define an orthogonal projector PfocI< such that 

-PFock^A-PFock = % F ° Cl \ (26) 

and describe the system in terms of "Ha = "Hp ock (g) and a supers election 
rule: a positive, normalized operator p represents a quantum state if and 
only if 

P = iWiW (27) 

Put differently, in this formulation of the EQDR theorem, the state carries 
the statistics. Positive, normalized operators that do not satisfy Eq. (27) do 
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not respect the fact that we are dealing with a system of indistinguishable 
particles, and must be discarded [32]. This reformulation of the kinematics 
of quantum fields/many-body problems makes the application of the EQDR 
theorem straightforward: everything works as in Section 3.1, whenever the 
state p satisfies the physical condition of Eq. (27). 

3.1.2. A basis- dependent version of the EQDR theorem 

The proof of the EQDR theorem did not rely on the use of any particular 
basis of Ha- Here we explore the consequences of a basis-dependent manip- 
ulation that we will later use in some examples, and helps to emphasize the 
connections to Hamiltonian lattice quantum field theory. 

As in Refs. [2, 3], we may consider arbitrary functions / of the boundary 
fields 4>o(x) and its conjugate momentum no( x ) — —i5/6<j>o(x). In such a 
case, / acts trivially on the bulk fields ip{x). Consider the orthonormal basis 
{|0o)} (without a tilde) to be an eigenbasis of / on the boundary F (meaning 
that / is a normal operator). The product states {|0o) <8> \4>) = 100^)}; i- e -> 
a direct product basis of states of V and the bulk A, form an orthonormal 
basis on which all physical states {\4>)} G Ha can be spanned, and are, 
trivially, eigenstates of / as well. For instance, for a situation in which 
/ = 7r (xi)TT (x2), we will choose the basis to be spanned by the set of 
product states, such that f\<fioip) = •^0(^1)^0(^2) |0o"0)- 

We need to evaluate Tr\(pf) = Tr^^^pf). A general matrix element 

(0o^|p/]0o^> =X) ^oV'bl^'X^'l/l^oV') = <^oV'bl0o^><0o^|/|0o^>,(28) 

as {\(po?p)} are eigenstates of /. Thus, 

(f) D = Tr {M} (pf) = Tr { ^ } (p diagona |/), (29) 
where the density operator Pdiagonai is defined as 

Pdiagonal = \MHM\p\MKM\, (so) 

{\m>)} 

such that ((potp I pdiagonai | <Poi ) ) > 0, reflecting the fact that only the diagonal 
elements of p matter in such a case. If we define a reduced density operator 
on A 

Pdiagonal [^] = ~T r {<j> } (Pdiagonal ) , (31) 
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then, clearly, 

Tr{^} (diagonal [^]) = T r {0o^} (Pdiagonal ) = ^{<f>o^}(p) = 1- ( 32 ) 

The average of / within the complete volume A (containing both the bound- 
ary and the bulk), 

(f)° = Tr w , } (Tr Wo} (p diagona |/)) = Tr w (p diagona |[<0] (/)$), (33) 

is the expectation value of / with fixed fields ip in the bulk A. In the last 
equality of Eq. (33) we used the fact that ^diagonal ['0] is invertible, 

(f)i = Pdiagonal['0] _1 Tr {0o} (p d i a g Ona |/). (34) 

We are interested in establishing bounds on (f) D . The latter may rep- 
resent an arbitrary correlation function on the <i-dimensional sub-volume T. 
Proceeding as in earlier subsections, we obtain 

(f)L n < iff < (f)t^ (35) 

where {^ max } and {VVn} represent the collection of fields in the bulk for 
which (f)^ is maximal/minimal. The expectation values (f)^ are weighted 
with a normalized positive weight Pdiagonai W\ • With the substitution of these 
fields, Hamiltonians such as -£^/, max [0o] and if,/, min [0o] (or corresponding ac- 
tions) will be functionals of only the boundary fields 4>q{x). Furthermore, 
similar to what we re-iterated earlier (and do so again now), these Hamilto- 
nians will be local in d dimensions (i.e., involve only local interactions) if the 
original .D-dimensional theory is local. As further noted earlier, this local- 
ity is to be contrasted with a brute force integration of the bulk fields x/i to 
obtain exact effective Hamiltonians i7 e ff[0o] within the sub- volume T, where 
arbitrary long-range interactions may be involved. 

We note that if we choose the quantity / in Eq. (35) to be the sym- 
metrized function / sym of Eq. (19) then the minimizing fields VWn (and 
similarly the maximizing fields t/i mM ) will be related to one another by the 
symmetry operations (d-GLSs) U.&. In particular, if ^ min (or -?/> m ax) transform 
as a singlet under the symmetry operations Ud then the averages (/sym)$ max min 
are computed with an effective local d-dimensional Hamiltonian (action) that 
respects the symmetry Ud- 
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3.2. The role of symmetries 

We next consider the consequences of symmetries on the bounds of Eq. 
(35) for systems for which the (correlation) functions /, the pertinent projec- 
tion operators, and symmetry generators that have their support on different 
spatial regions commute with one another (and thus Eqs. (37, 39) below will 
be satisfied). We note that with the fixed bulk fields Vw x and ^ min that 
form the bounds in Eq. (35), the d- dimensional Hamiltonians -£^/> max [0o] and 
H^ min [(f)o] (or actions) will admit a global d-GLS symmetry on the region T 
if the original D- dimensional theory exhibits a d-GLS on that region. (This 
is so as by virtue of the existence of d-GLSs, the bulk fields ip must act 
as external non-symmetry breaking fields in -£^/, max [0o] and if^ min [0o] or any 
other Hamiltonian formed by fixing i/j: if^[0o]-) The inequalities in Eq. (35), 
leading to the concept of dimensional reduction in the weak sense of upper 
and lower bounds, are quite general for systems displaying d-GLSs. As noted 
earlier, the operators / are quite arbitrary and may be chosen to be arbitrary 
n-point correlation functions. When {a^}™ =1 G T, the exact n-point corre- 
lators within the D- dimensional theory are bounded, both from above 
and from below, by correlation functions in local theories in d dimensions 
(theories defined on T alone) that exhibit d-GLSs, 

G\l r , d ({xi}) < gSL ; d({^}) < G [ :J pend ({x t }). (36) 

Equation (36) fleshes out a corollary of the EQDR theorem for systems with 
d-GLSs. Similar to our discussion of Eq. (19), we may choose the correlators 
G^ to be symmetrized relative to the symmetry generators. This may ef- 
fectively set the averages ((j^) U p P er/iower;d to be those computed with a local 
d-dimensional theory that respects the symmetry generators at hand. 

It is important to emphasize the reason for being able to freeze out the 
bulk degrees of freedom in obtaining the inequality (35), while maintaining 
d-GLSs in the bounding theories if the original theory exhibits d-GLSs. Let 
us define the projection operator which has its support only on degrees 
of freedom in the bulk A, while the lower-dimensional symmetries {U d } have 
their support on fields in the sub-volume V. Therefore, 

[P^,U d ]=0 (37) 

as A PI r = {0}. As U d are symmetries of p, we then have that 

[PtpPtM = P^pMPi, = 0. (38) 
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The expectation value of Eq. (34) is that of f(4>o) with the density oper- 
ator p computed over the projected density operator p^ = P^pP^ when it is 
invariant under Ud- This is so since 

[/(0o),P^] = O (39) 

(as, once again, the two operators / and P^ have their support on T and A, 
respectively). This commutation relation implies that, e.g., in the canonical 
ensemble, 

Uh Tr^y^e-WP^) Tr W} (e~^) " 1 ] 

It is interesting to highlight those cases in which / is not invariant under 
low- dimensional GLSs (d = 0, 1 discrete symmetries or d = 0, 1, 2 continuous 
symmetries) [2] in a system with interactions of finite range and strength. 
Then, as such symmetries cannot be broken in the low d-dimensional sys- 
tem, the expectation value of / must vanish at finite temperatures [2, 3]. 
The case of d = corresponds to Elitzur's theorem for usual gauge theories 
(with local (d = 0) symmetries) [33]. One-dimensional systems with inter- 
actions of finite range and strength do not exhibit spontaneous symmetry 
breaking at non-zero temperatures. Similarly, two-dimensional systems do 
not exhibit spontaneous symmetry breaking of continuous symmetries. Thus, 
Elitzur's theorem can be extended to such cases [2, 3]. Similarly, if a low d- 
dimensional correlation function decays (exponentially or algebraically) with 
distance then by virtue of our inequalities the related correlation function in 
the D-dimensional theory is similarly bounded. The physical engine behind 
the absence of symmetry breaking (as well as more exotic high- dimensional 
fractionalization) in high-dimensional systems with these symmetries is that 
the energy-entropy balance in these D-dimensional systems is essentially that 
of d < D dimensional systems [2, 3, 34]. 

Our EQDR theorem and the ensuing consequences of d-GLSs may seem 
rather formal. To elucidate the physical consequences of these symmetries, 
we turn to the example of the square lattice compass model of Eq. (2) [2, 3]. 
Let us first consider the consequences of the EQDR theorem of Eq. (35) 
for this system. If we consider the boundary T to be a line (i.e., a d — 1 
dimensional object) along a lattice direction (say, the e.\ direction) in the full 
two dimensional (D = 2) system of i^compass; then the correlation function 
(f) D = (&x a x') between two sites x and x' that lie on the line T is bounded 
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(both from above and from below) by the same correlation function in one- 
dimensional systems. These one-dimensional systems in which we compute 
(/)$ max and (/}$ min may be seen to have an effective Hamiltonian that is 
composed of nearest-neighbor Ising type interactions along T augmented by 
external transverse fields (originating from the a = y contributions in Eq. (2) 
once we fix the bulk fields (spins in this case) on all sites y G A, i.e., outside 
the line T, to a set of particular values {^w} or {ip max } (or {cr^. min }, {cr^; max } 
in this case). That is, the averages, as in Eqs. (16, 17), are performed for 
d = 1 systems. For each vertical or horizontal line T we have the symmetries 
of Eq. (3). On any given d = 1 line T, these symmetries will appear as global 
(d = 1) symmetries. The external (transverse) fields ip do not spoil these 
symmetries. 

As in any such transverse field type Ising chain, there is no long-range 
order at any finite temperature and coupling, the correlation function of 
(< 7 x a x') tends to zero as sites x and x' become well separated. As noted 
in Eq. (3), H cornpass is invariant under the following d = 1 (Ising type) 
symmetry operations of U& = Ylxet a x where Y denotes any (d = 1) line 
parallel to the ex-axis . (A similar set of symmetries applies for lines parallel 
to the e2-axis with the x and y Pauli matrices interchanged.) In this case, 
we may, in particular, choose V to coincide with the region V where the 
symmetry operates. On the line T, we have (as stated earlier) in the bounding 
theories on both the left and right-hand sides of Eq. (35), a short-range 
one- dimensional Ising model compounded by transverse fields which do not 
break the d = 1 Ising symmetry above. This absence of symmetry breaking 
constitutes an analogue and extension [2, 3] of Elitzur's theorem that pertains 
to local gauge (i.e., d = 0) symmetries [33]. When Ji ^ J 2 , this analysis 
implies that the D = 2 compass model does not display long-range order 
at any finite temperature. At the isotropic Ji = J 2 point an additional 
d = D = 2 discrete symmetry appears which can be broken. 

4. Exact dimensional reduction: general scope 

In the next two sections, we discuss two related forms of exact dimensional 
reduction. The two cases studied are those that appear in (a) various limits 
such as large n in 0(n) vector theories or high temperature (or weak coupling) 
and (b) dimensional reductions via exact dualities that do not appear in 
some limit of the theory but can rather be established by considering the 
relevant degrees of freedom (so called bonds) that define the Hamiltonian (or 
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action) of the theory [1, 18]. Dimensional reductions of type (a) above can 
be established by considering the density of states (DOS) as we will explain 
below. Systems such as those discussed in Section 6 and Kitaev's and Wen's 
models [19, 35] illustrate dimensional reductions of type (b) [2]. In such cases, 
similar to [1, 18] a (unitary) transformation exactly maps the two systems 
onto one another and preserves their spectrum. 

5. Exact dimensional reduction from the density of states 

We will now illustrate how a technique based on the DOS as applied to 
n-component vector theories (i.e., "O(n) theories") will enable us to estab- 
lish dimensional reductions in various instances. This includes some limits 
of 0(n) vector theories such as that of large-n or that of weak coupling or 
high temperatures. In a Feynman diagrammatic framework, these relate to 
situations in which we only consider tree level diagrams. As emphasized ear- 
lier, although symmetries lead to bounds on correlators such as those of Eqs. 
(35, 36), they generally do not lead to an exact dimensional reduction per se. 
That is, the thermodynamic functions in systems with a d-dimensional sym- 
metry are, generally, not those of canonical d-dimensional systems but rather 
reflect the full dimensionality of space-time. There are, however, very gen- 
eral situations in which the dimensional reduction is guaranteed to be exact 
thermodynamically (i.e., the thermodynamic functions of the D-dimensional 
systems are identical to those of d- dimensional systems). 

We next demonstrate that in the large-n limit of all vector theories, the 
spatial dimensionality of the system can be changed at will if the total DOS 
profile is kept the same. 

5.1. Preliminaries and statement of the DOS theorem 

We now focus on lattice field theory Hamiltonians of the form 

x,y x 

where (i) the pertinent vector fields (0 = (0i, <f> 2 , • ■ • , 4>n)) in Eq. (41) have 
n ^> 1 (large n) components and (ii) the strength of the quartic coupling in 
Eq. (41) is set by u = 1/n. In what follows, we denote the Fourier transform 
of the translationally invariant J(x — y) by J(k). 
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Given the general form of Eq. (41), we define the DOS via 

0o(e) = J ^0j^o(e,k), (42) 

where the local DOS is g (e, k) = 5(e — J(k)). In Eq. (42), the fc-space inte- 
gration is performed over the first Brillouin zone for lattice theories. In the 
regularized continuum field theoretic rendition of Eq. (41), the integration is 
performed over all real vectors k. As we will briefly discuss, the DOS of Eq. 
(42) is that of the non-interacting variant of Eq. (41), i.e., one with u — 0. 
With these preliminaries, we are now ready to state our next theorem: 

Theorem 2 (DOS): The large-n vector system of Eq. (J,.!) in D dimensions 
and another large-n vector system of the same form in which the real space 
interaction kernel is replaced by another kernel that resides in D' dimensions, 
J(x — y) — > J(x' — y'), will have identical free energies if their DOS are equal, 
i.e., if 

J d D k 5(e - J(k)) = J d D 'k' 5(e - J (Is!)). (43) 



The proof of this assertion follows from rather standard results that we 
show below. We first remark that this constitutes an extension of the trivial 
variant of this statement for the non- interacting system of u = 0. In such a 
non-interacting system, Eq. (42) provides the exact DOS for a single mode 
(and, by extension, to the collection of decoupled modes that constitute 
the non-interacting system). Identical forms appear for other systems of 
decoupled particles as in, e.g., the DOS for free electrons wherein the kernel 
J(k) is replaced by a band dispersion e(k). For the lattice systems that we 
consider, each Fourier mode constitutes such a free particle; the number of 
such particles (which we generally generally denote by N) is equal, for the 
lattice systems that we focus on, to N s (= |A|), the number of sites of the 
lattice. The total DOS is given by 

g{E) = J de x ■ ■ ■ J de N 5{E - e x - e 2 e*). (44) 

Two systems are dual to each other and share the same spectrum if and only 
if their densities of states g(E) are identical [1]. Thus, in the non-interacting 
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limit, when Eq. (43) holds, the two corresponding systems are trivially dual 
to one another. The special character of the large-n vector theories renders 
this result also precise also when interactions are introduced (u—l/n). 

The crux of our proof of the DOS theorem for large-n vector theories, and 
more general possible extensions thereof to other arenas, is that the partition 
function depends solely on the DOS of the system 

Z = Z[g(E)]. (45) 

Indeed, by construction, as Z = Tr\(exp(—(3H)) = f dE g(E) exp(— j3E), 
the partition function is the Laplace transform of the full exact DOS of the 
system g(E). As such, any transformation that leaves the density of states 
of the full system g(E) invariant automatically ensures that the free energy 
is the same. Two different systems of free particles (i.e., quadratic systems) 
with an identical single particle DOS go{e) (that determines the full DOS 
of the system) must, by Eq. (45), display an identical density of states 
g(E) for the entire system. In such (free-particle type) Gaussian systems, 
the partition functions can be further easily computed in the presence of 
external sources. That leads to a quadratic coupling between the sources. 
In interacting systems, finding such a description with an effective density of 
states <7eff(e) is, generally, far more complicated. The dependence of the full 
interacting system on a single density function is reminiscent of the density 
functional theory that has been of immense success in solving quantum many- 
body systems. The Hohenberg-Kohn [36] theorems that encapsulate this 
approach illustrate how the system properties may be effectively captured by 
a single density function. A physically appealing property of the dimensional 
reduction that we outline here is that, while preserving the DOS functional 
(set by go(e) in large-n vector systems), we may map a local system (one with 
interactions that decay with distance) in some spatial dimension onto another 
local system residing in a lower spatial dimension. As outlined above (and 
below), this preservation of the spectrum, i.e., the DOS, while maintaining 
the local character of the theory by performing such a map constitutes the 
key physical underpinning of dualities [1]. 

We will next review some results concerning the large-n vector theories 
to explicitly illustrate how the result of the DOS theorem is manifest. 

5.2. Free energy of large-n vector systems 

The DOS theorem follows nearly immediately from a diagrammatic cal- 
culation [37] as well as the well known equivalence between large-n vector 
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theories to the spherical model [38]. 

Specifically, the free energy of (the to be reviewed) Eq. (53) wherein fi is 
implicitly determined by Eqs. (51, 52) is manifestly dependent only on the 
DOS go{e). This result similarly follows from diagrammatic considerations. 
Physically, as will elaborate later on in an example, the quintessential low- 
energy features of the theory are determined by the soft modes. 

We review several known key features. In the large-n limit, "hard spin" 
models (those with u ^> 1 in Eq. (41) in which the fields have a fixed hard 
constraint on their norm, i.e., \(j)\ = 1) tend to behave as their "soft spin" 
counterparts (with u ~ 1/n). Furthermore, the large-n rendition of Eq. (41) 
is, in many respects, equivalent [39] to that of the spherical model [40]. In 
particular, the Helmholtz free energy density F n associated with the 0(n) 
system of Eq. (41) is, in the limit n — > oo, identical to the free energy density 
-Spherical of the spherical model version of the 0(n) theory, i.e., F n = F sp h e ricai- 
On a finite size lattice, the n-component spherical model [41] is defined as a 
theory of 0(n) fields with Hamiltonian 

H = I E J ( x - y)^ x ) ■ ttv) = 2F E wiwr, (46) 

x,y S k 

augmented by the spherical constraint, 

X 

The sum over wavevectors k is constrained to those in the first Brillouin zone. 
In the thermodynamic limit, this sum may be replaced by an integral. We 
may implement the constraint of Eq. (47) via a Lagrange multiplier (u/2). 
Thus, in the large-n limit, the quartic term of Eq. (41) may be replaced 
by the spherical constraint of Eq. (47) and the theory is defined by the 
Hamiltonian 

^spherical = \ ^ J(x - y)<f>(x) ■ <f>{y) + | ^ 4>\x). (48) 
x,y x 

In order to compute averages, we augment //spherical by local fields h(x) 
that couple linearly to the vector field 4>(x) and take functional derivatives 
of the resulting generating functional Z relative to these fields. Thus, in 
the spherical model, or lowest (zeroth) order 1/n computation, a Gaussian 
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integration over the fields leads to the corresponding canonical generating 
functional 



^spherical = D<f) exp(- f3[H spberica \ h(x) ■ <f)(x)]) 

x 

rj r 2ti . p ^h(k)-h(-k). 

= n ^ "pi-, E ^tttI. m 

where D(f> is a shorthand for the lattice measure W^. 3 d(j)(x). By differentia- 
tion relative to the fields h(k) (or, alternatively, by the use of the equiparti- 
tion theorem), it is seen that 

<w k)|2 > = mkrjy (50) 

Inserting Eq. (50) into Eq. (47) and accounting for a possible condensate 
(detailed below), we arrive at the self-consistency equation 

= Jw?WnT, + m '- (51) 

The last term represents, at temperatures below T c (if a critical temperature 
exists), the condensate of modes q that lie on the manifold of minimizing 
modes M. That is, J(q 6 M) = min k {J(k)} and in Eq. (51), 

As is evident from Eq. (49), a permutation P of the wavevectors k (i.e., 
a transformation J(k) — >■ J(?k) = J(k')) will leave the product over all k 
invariant. The generating functional depends only on the density of states 
g(e). Then, it is clear that the zero field (h = 0) free energy is given by 

Spherical {h = 0) = -/T 1 In Z spheric ai (h = 0) = ^ de g {e) ln[/3(e + //)], (53) 

which clearly demonstrates the DOS theorem of Eq. (43). 

These results also follow from the direct diagrammatic expansion of the 
theory of Eq. (41) to lowest order in (1/n) (i.e., to order 0((l/n)°)) [37, 42]. 



23 



5. 3. Examples of dimensional reduction in large-n theories 

We should stress that the dimensional reduction that arises in the large-n 
theory is general and not limited to these examples. Below, we will review 
facts pertaining to the DOS in canonical systems with a quadratic disper- 
sion about the minimum of the energy in /c-space and employ these in two 
examples of dimensional reductions that highlight quintessential aspects of 
the DOS theorem. Insofar as the simple theory of Eq. (2), an explicit 
calculation in the large-n limit illustrates that this system is identical to a 
one- dimensional system in the same limit. We now provide two simple (yet 
slightly less trivial) examples. 

Example 1: Consider the case of nearest-neighbor interactions on a 
hypercubic lattice in D dimensions in which the lattice constant a is set to 
one. In such ferromagnetic lattice large n model is given by Eqs. 

(46, 48) with J(x — y) = — J if sites x and y are nearest neighbors and is 
zero otherwise. In Fourier space, the corresponding kernel is given by 

D 

J(k) = 2 J ^(1 - cos ki) = JA(k). (54) 
e=i 

As seen, here, A(k) denotes the lattice Laplacian in D dimensions. In the 
low energy or the small k (long-wavelength) limit A(k) — > k 2 . The number 
of points in /c-space that have a given value of J(k) = e is set by the surface 
area of a sphere whose radius is given by (e/J) 1 / 2 . Thus, the DOS scales as 
go(e) ~ g( £) - 1 )/ 2 . In Eq. (54), the quadratic minimum (for J > 0) occurs at 
a single point in /c-space, k = 0. 

Let us start by considering the case in which the kernel J(k) of Eq. (41) 
has its minimum on a manifold M of dimension <i min and that J(k) exhibits 
a quadratic dispersion about this minimum. In such cases, the DOS is that 
of a D' — (D — drain) dimensional system. That is, to leading order 

9o(e) ~ e^'^l\ (55) 

The dependence of the generating functional solely on the density of states 
go(e) allows all translationally invariant systems with a particular momentum 
space dispersion to be mapped onto a one-dimensional system. Towards 
this end, we define a kernel J e ff{k) of an effectively equivalent one- dimensional 
system (i.e., that with a scalar k) by 



/ 
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The above relation ensures that the DOS and consequently the generat- 
ing functional are preserved. That is, the generating functional of the D- 
dimensional system with the kernel J(k) is identical to that of a one-dimensional 
system with the kernel J e ff(k). 

As an example, we now consider the nearest-neighbor system of Eq. (54) 
in two dimensions (D = 2). We will now map this system onto a d = 1 
dimensional system. Towards this end, we write 



dJ e f( 


-l 


r dk x f 


dk 


J eff (k)=e 





t jfTwuhw^ = 9o(£) ' (57) 



with e = e/(2J). The second equality in Eq. (57) trivially follows from, 
e.g., the substitution v = cos&i (after w = cos/c 2 is integrated over and 
consequently set by the delta function). Thus, 

fc(Jeff) = / deg {e). (58) 



Equation (58) may be inverted and Fourier transformed to find the effective 
one-dimensional real space kernel J e ff{x — y). By inserting this back into Eq. 
(41), we complete the mapping of the two-dimensional nearest-neighbor sys- 
tem onto a one-dimensional system. In a similar fashion, within the spherical 
(or equivalently the 0(n — > oo)) limit, all high-dimensional problems may 
be mapped onto a translationally invariant one-dimensional problem. It is 
important to re-iterate, in light of the example provided above, that the 
lower (one-) dimensional model albeit having longer range interactions than 
its higher dimensional counterpart (in fact, scaling as V(r) ~ r~ 2 for large 
r = \x — y\), the resulting theory is nevertheless local and physical. 

Example 2: As another more natural example where dimensional re- 
duction occurs, we consider next a situation wherein J(k) attains its global 
minimum on a spherical shell |k| = q > 0. When J(k) attains its minimum 
on a spherical shell, then in an expansion of about any such point, there 
is only one direction in fc-space (the radial direction) along which the dis- 
persion may be quadratic and (D — 1) tangential directions with soft modes. 
The DOS is determined solely by the massive radial direction, and this sys- 
tem behaves as a one-dimensional one. The correlation functions are given 
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by derivatives of the generating functional relative to h in which the limit 
h — y is taken. If the low energy dispersion is given by Eq. (55), then the 
free energy (in which only states in the vicinity of the minimizing manifold 
are considered as is appropriate at low temperatures) will have the canonical 
form of a D' = D — d m i n dimensional system. That is so as in the coordinate 
system locally tangential to the minimizing manifold in /c-space, there are 
D directions on which the dispersion J(k) has a quadratic dependence and 
c? m i n directions on which it has no dependence. Thus, in the large-??, limit, 
not only are some correlation functions of a form characteristic to that of a 
(lower) D'-dimensional system as we will show but also the thermodynamic 
functions themselves are also characteristic of a low- dimensional system. Us- 
ing Wick's theorem for the quadratic large-?? theory, the correlation functions 
between 2m fields are given by 



Cast in terms of symmetries, the operation that enables us to reduce 
the D- dimensional problem onto a /^'-dimensional one is the presence of a 
non-trivial permutational d m \ n - dimensional symmetry of the spectra. That 
is, any permutation P (whose minimal non-empty set spans manifolds M of 
dimension <i m i n ) that maps wavevectors k with the same energy (J(k)) onto 
each other is a symmetry. The full <i m i n -dimensional symmetry enables us to 
reorganize the labeling of states in /c-space such that the resultant problem 
after a permutation corresponds to a Z)'-dimensional one (or more precisely, 
at low temperatures \M\ copies of a D'-dimensional system where \M\ is the 
number of minimizing modes). We briefly elaborate on this. We can choose 
permutations contort the general G? min -dimensional equipotentials of constant 
J onto a set of flat <i m i n - dimensional hyperplanes. These represent a stack of 
decoupled .D'-dimensional systems (the stacking occurs along d min directions 
along which there are no interactions: there is no mode dispersion). 

5.4- Weak coupling or high temperature 

We now discuss an extension of the DOS theorem of Eq. (43) to an- 
other class of systems. As illustrated in [43], in the high temperature (or 
weak coupling) limit of any 0(n) theory (with an arbitrary finite number of 
components n) of the form of Eq. (41), the free energy is given by 
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(59) 




(60) 
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Thus, similar to the case of the large-n vector theories, all that matters in Eq. 
(60) is the DOS g{E). In particular, any system can have its dimensionality 
reduced. There are numerous systems that, at the high-temperature (small 
j3) limit have an identical free energy. A subset of these systems become 
lower dimensional both (i) in the limit of small ft (for finite n) and (ii) the 
limit of large n for finite j3. For instance, a one-dimensional Ising chain 
with a momentum space interaction kernel given by Eq. (58) with an Ising 
version of Eq. (46) is identical, in the limit of high temperatures, to the 
two-dimensional nearest-neighbor Ising model. 

5. 5. Dimensional reductions and density of states functional theories 

Beyond the particular confines of the large-n vector theories and high 
temperature/weak coupling limits, the likes of the DOS theorem might hold 
elsewhere. We now make general remarks about the general premise and 
viable extensions of the DOS theorem. We first couch our specific consid- 
erations thus far within a simple standard framework. As highlighted in 
Eq. (45), in these theories, a single (DOS) functional determines the system 
properties as the DOS determines the generating (or partition) function Z. 
As it so happens in the 0(n) theories that we looked at the free energies at 
hand were a DOS functional set by the interaction kernel in /c-space. The 
DOS theorem would hold without any change for more general functionals. 
We remark that density functional theories with correlations set by the local 
DOS in real space representations (au lieu of the fc-space dependence that 
investigated in the large-n theories) have been advanced in solid state sys- 
tems [44]. Just as in many-body systems, it well may be hard to determine 
the DOS functional apart from easier cases such as that of large n. 

We speculate that a dependence on the metric g^ix) in gravity the- 
ories might play an analogous role to the dependence on the interaction 
kernels that we investigated above for large-n vector theories, and that ef- 
fective metric functional theories might exist by an extension of standard 
methods of density functional theories [36]. In such instances, computable 
quantities would be set by functionals of the metric. We would like to em- 
phasize, though, that an equivalent of the Hohenberg-Kohn theorems should 
be proved first. Had this proof be non-existent then there is no possibility 
of an exact duality (i.e., an isomorphism) between quantum many-body or 
non-relativistic field theories and semiclassical gravity theories. It could still 
happen that a duality emerges in some appropriate limit. 
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We now step back to highlight the specifics of the core features of the 
derivation of the DOS theorem and the applications that we provided. Our 
treatment was centered on the Gaussian type character of large-n theories. 
In Eq. (49), we introduced sources. The functional derivatives relative to the 
source fields h yield the correlation functions. Written in a more general form, 
Gaussian integration leads to a free (i.e., quadratic order) form in the source 

(2) 

fields which are coupled via the correlation function G . The standard 
large-n results that we invoked and expanded on may be more generally 
linked to a standard perturbative framework. Within a general perturbative 
framework the kernel [J(x — y) + /i] of Eq. (48) coupling the fields <p(x) 
and 4>(y) may, more generally, be replaced by a zeroth order (i.e., quadratic) 
two-point kernel [(Gq 2 '') -1 ]. In such the generating functional is 



Z (/3,/i) = Z (/3,0)exp 



~ / d D xd D y h(x)(PG { *\x,y))h(y) 



(61) 



As is well known, for a general Hamiltonian (or action) density, which 
is the sum of a quadratic and an interacting part "Hint, the corresponding 
generating functional becomes 



2(/3, h) = exp 



J d D x f3 H mt {-6/6{f3h{x))) 



Zo(P,h). (62) 



Thus, whenever dimensional reduction is exact on the quadratic level, it may, 
in certain cases (e.g., that of Example 2 above) enable effective dimensional 
reductions when a perturbative expansion about the quadratic theory is valid. 

There is, generally, no expression for the lowest-order partition function in 
the presence of sources, e.g., Eq. (49), in terms of the single particle DOS. In 
the presence of sources, the partition function depends on the full dispersion 
J(k) in fc-space and not solely on the single particle DOS go{e) (nor the 
full DOS of the entire system g(E) in the absence of sources). This feature 
underscores the larger information content that may be generally present in 
high-dimensional fc-space when sources are present, vis a vis that embodied 
in the system dependence only on the single energy coordinate e (or E). We 
employed this redundancy in the higher-dimensional /c-space in order to freely 
move from a system in D dimensions to a system in a lower dimension when 
sources were not present (Example 1). In situations in which two dispersions 
J(k) and J(k') in different dimensions (say D and D', with D > D') exhibit 
identical behavior in the vicinity of their minima then the partition functions 
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in the presence of sources and (the corresponding correlation functions) may 
be identical. This path was realized in Example 2 above. In such a situation, 
the high-dimensional function J(k) is constant, i.e., does not disperse, for 
variations in k along D — D' = d m - m directions, while along D' components of 
k the function J(k) is identical to J(k'). As J(k) captures, to lowest order, 
the correlation function G , two such systems exhibit identical correlations. 
Moreover, the identity of the dispersions J(k) and </(k') ensures that the 
corresponding DOS are the same and, thus, the free energies are identical. 

We now briefly remark and speculate on possible connections between the 
above results and those suggested by AdS-CFT correspondences [10]. The 
relation of Eq. (43) enables a construction of candidate gravity duals for 
0{n) theories. Recently [45], an attempt has been made to find the gravity 
dual of an 0(n) vector theory using a renormalization group scheme. It was 
found that an 0(n) vector theory corresponds, via an AdS-CFT type corre- 
spondence, to gravity in the Arnowitt-Deser-Misner (ADM) representation 
[46] with two sets of conjugate fields. Perusing the results of Ref. [45], we 
note that in the large-n limit, the dual ADM theory is quadratic (as it must 
be) just as the original large 0(n) field theory to which it is dual. With the 
results of that we derived in the current section, it may be straightforward to 
check whether the relation of Eq. (43) is satisfied in this case. This is so as 
our relation of Eq. (43) enables the appearance of dual equivalent theories 
of many sorts, not only dual gravity theories to standard CFTs, that exist 
in one higher dimension. 

The ground state entropy of vector theories in the large-n limit has a 
holographic character, i.e., it scales as the surface area of the minimizing 
manifold M in fc-space 

S ~ (63) 

As in our discussion of Example 2, in Eq. (63), <i m i n is the dimension of the 
surface M of minimizing modes in /c-space [47] . The proof of this assertion 
is trivial. Any state 0(x) = J^qeM ^(^0 6xp[zq • x] (with <f>(— q) = 0*(q) to 
ensure that 0(x) is real) whose sole non- vanishing Fourier modes are those 
for which J(k) attains its minimum is a ground state. The manifold M 
has dimension d m \ n and thus there are 0(L dmin ) minimizing modes q and 
corresponding Fourier amplitudes 0(q). These amplitudes have to adhere to 
the single global normalization constraint of Eq. (47). This leaves a total 
of 0(L dinin ) free amplitudes {0(q)} and to the sub-extensive holographic-like 
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entropy of Eq. (63). For finite-n systems, not only the density of eigenmodes 
in fc-space is important and the above correspondence breaks down. 

6. Exact dimensional reduction from bond algebras 

In this section we will rely on a new bond-algebraic theory of dualities 
[1, 18] to study exact dualities that connect models of dimension D to models 
of dimension d < D. These dualities afford especially transparent examples 
of exact dimensional reduction. Dualities have been the subject of intense 
research efforts for some years now since the AdS-CFT correspondence is 
conjectured to be a particular example of a large class of dualities that ef- 
fect a change of dimension [13, 14]. The subjects of TQO and fault-tolerant 
quantum computation are also rich in models in D > 2 dimensions that are 
dual to d = 1-dimensional models. This is not surprising if one notices that 
these .D-dimensional models typically display d= 1-GLSs [2] and, as we have 
seen already, d-GLSs are closely linked to effective dimensional reduction. 
However, the precise connection between d-GLSs and exact dimensional re- 
duction requires further discussion. We would like to understand whether 
1) models that display d-GLSs are always exactly dual to lower-dimensional 
models, and conversely, whether 2) D-dimensional models dual to d < D- 
dimensional ones must display d-GLSs. That d-GLSs do not always imply 
exact dimensional reduction will be shown by a counter-example in the next 
section. As for 2), we do not have a full understanding yet. 

There is a different aspect to the connection between effective and exact 
dimensional reduction that also calls for clarification. The bounds derived 
in Eq. (36) are not saturated in general, meaning that the upper and lower 
bounds G^pp erd and G\^ erd may be different. They nevertheless afford a 
sharp characterization of effective dimensional reduction, and from this per- 
spective, exact dimensional reduction may be understood as the special case 
in which inequalities of the form of Eq. (36) become equalities. It is natural 
to conjecture that this happens if and only if there is an exact duality that 
maps the D-dimensional model to an explicitly lower dimensional model. 
The reason is that a duality amounts to a unitary equivalence [1] . It follows 
that the correlators of several appropriately defined quantities in a model 
of dimension D > d dual to a model of lower dimension d must show ex- 
act d-dimensional behavior, because they can be mapped to and represented 
by corresponding duals in the low (d) dimensional system. This is a fairly 
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simple proof of the "if part of the conjecture. At present, we do not know 
beyond some specific examples whether the "only if part holds as well. 

Regardless of their connection to effective dimensional reduction, exact 
dualities that change dimension will be studied here also for their own sake. 
We will see that they represent a delicate phenomenon that results from a 
very peculiar, hard to characterize structure of interactions (bond algebra in 
the language of [1, 18]), and that dual models of different dimension can only 
be modified in limited ways without destroying the duality that connects 
them. We think that these patterns visible in concrete examples should be 
kept in mind when gauging the merit of conjectured dualities that change 
dimension, and most importantly, when trying to extend these conjectures 
in non-trivial ways. 

We will focus both for simplicity and for their own merit on lattice models 
of TQO. If a system displays TQO, its state cannot be solely characterized 
by local measurements. Rather measurements of topological quantities are 
needed too [20, 19]. Systems characterized by topological (non-local) order 
are, in principle, robust to local perturbations and thus able to store quantum 
information for long periods of time. Such systems would be ideal candidates 
for realistic quantum memories for quantum information processing and com- 
putation. Unfortunately, many early models harboring TQO were shown [2] 
to be dual to one-dimensional Ising chains with short-range interactions, a 
result that suggests that the autocorrelation times in these particular models 
are system-size independent. Moreover, at equilibrium all initial information 
becomes random and is lost at any non-zero temperature, a phenomenon 
dubbed thermal fragility [2]. Here we will show that this property is also 
shared by several new models that have been devised and investigated since 
the early days of TQO. They are: 

1. a D = 2 quantum model on a honeycomb lattice generalizing the one 
introduced in Ref. [48], 

2. D = 2 topological color codes [49], and 

3. a D = 3 model on an fee lattice introduced in Ref. [50] and further 
investigated in Ref. [51]. 

These models are all dual to one-dimensional Ising chains or ladders with 
appropriate boundary conditions, and appropriately distributed transverse 
fields. They display exact dimensional reduction, and, as a consequence, 
memory times that are system-size independent at all temperatures. 
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In closing, let us notice that the techniques and ideas to be discussed next 
are perfectly suitable [1] for studying models of lattice Hamiltonian quantum 
field theory [52]. This observation may become specially relevant if and when 
lattice renditions of the AdS-CFT and related conjectures become available 
for direct numerical testing. 

6.1. Extended toric code on a honeycomb lattice 




Figure 1: A honeycomb lattice described as a triangular Bravais lattice with a basis. 

The model that we are going to study in this section is a variation of 
a model introduced in Ref. [48]. It features 5=1/2 spins on each site of 
a honeycomb lattice, and shows exact dimensional reduction despite having 
non-commuting operators in its Hamiltonian. 

The honeycomb lattice can be described as a triangular Bravais lattice, 
generated by the primitive vectors a x = e x and a 2 = (e± + \^3e 2 )/2 (e% and 
e 2 are the standard unit vectors of a square lattice), together with a basis 
described by the two vectors and b = [e.\ + (e 2 /v / 3))/2, see Fig. 1. The 
sites of the lattice are spanned by the sets of points m = m 1 a 1 + m 2 a 2 , 
and the set m' = m + b, with m 1 ,?^ 2 £ Z. Imposing periodic boundary 
conditions amounts to identifying two sites rai,m.2 to be the same if their 
coordinates m 1 ^ are equal up to a multiple of the period L M , \i — 1, 2; and 
we identify two sites m^, m' 2 if m^ — b, m' 2 — b satisfy the same condition. The 
result is a honeycombl lattice A^ ex with the topology of a two-dimensional 
torus. Since we can associate each elementary hexagon to a site m in the 
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triangular sublattice, we can use m as an index to label the two plaquette 
operators 

tt t z x y z x y 

yy{l,m) = <J m <T m' a rn+a2 a rn'+a2-a 1 a m+a 2 -a 1 a m'-a 1 ' 

W(2,m) = ^m^^m+aj^m'+az-a^m+^-a^m'-a! • (64) 

where with a — x,y,z are Pauli matrices. Notice that any two of these 
plaquette operators commute, [W( MiTn ), W(„ )n )] = 0. 

The honeycomb extended toric code (HETC) model that we will study is 
specified by the Hamiltonian 

Li-l L 2 

#HETC = - Yl Yl J V W (»,rn) ~ k °™ ( 65 ) 

//=1,2 rnGAF 771*=! m 2 =l 

hex 

illustrated in Fig. 2 for a lattice with periods L\ = 3, L 2 = 2. The Hamil- 
tonian acts on a Hilbert space of dimension (2^) 2 = 2 2Ns , with N s = L\L 2 . 
If h = 0, i^HETC reduces to the model introduced in Ref. [48]. Notice that 
the magnetic field h is staggered, meaning that it is zero on all sites m', 
and that it is mostly homogeneous, except for the (one-dimensional) circle of 
sites Lidx + m 2 a 2 , m? = 1, • • • , L 2 , where it is also zero. The Hamiltonian 
-^hetc is self-dual under the exchange J-y <H- J 2 , but we will not show this 
directly, as it will become clear after the duality to a d = 1 model. 

We will use the new approach to dualities developed in Refs. [1] to look 
for dual representations of -£/heto Consider the elementary operators W^ m ) 
and <j^j, and only those operators. It is convenient to have a uniform denom- 
ination for them, and so we agree to call them generically bonds [18]. The 
bonds of a Hamiltonian (like -Hhetc) generate a model-specific bond algebra 
of interactions that has typically a very different structure as compared to 
the algebra of elementary degrees of freedom (spins S = 1/2 in the case at 
hand). So it was proposed in [1] to look for dualities as structure-preserving 
mappings (homomorphisms, or isomorphisms in the absence of gauge sym- 
metries) of bond algebras local in the bonds. To find a dual representation of 
-^hetc we need to characterize its bond algebra, and find a model -^hejc ^ na ^ 
looks different but has in fact a bond algebra isomorphic to that of -Hheto 

Since each bond a z anticommutes with four plaquettes, -Hhetc's bond 
algebra is non-commutative. The bonds commute with each other, 

and satisfy the constraints 

]JW M = t, /i = l,2. (66) 
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Figure 2: Honeycomb lattice with L\ = 3, L-2 = 2. and periodic boundary conditions 
(twelve distinct sites). The top panel shows the plaquette bonds VF(i iTrl ) introduced in 
Eq. (64), numbered from one to six, and the four transverse fields a z as green crosses 
(the spin labelled as three is not part of the Hamiltonian, but it appears in one of its d- 
GLSs). The bottom panel shows the plaquette bonds W^ m )- The green (blue, red) dots 
stand for spins a z (a x , a y ). The middle panel represents the dual, one-dimensional ladder 
model defined in Eq. (67). The duality mapping defined in Eq. (68) is indicated by 
corresponding labels in the illustrations of both models. For example, the plaquette bond 
in the bottom panel labelled as 1 maps to the link ct| jcrf 2 indicated by a green horizontal 
bar, and labelled by 1 as well. The dual model satisfies periodic boundary conditions. 



Keeping these relations and constraints in mind, one can check that the bonds 
for a one-dimensional ladder of spins with periodic boundary conditions sat- 
isfy identical relations and constraints, as illustrated in Fig. 2. It follows 
that -£?hetc is dual to 

N s L 2 -2Li-l 

J V a l,,m a l,,m+l - h ^2 CT l,i+sC£*-l) a 2,i+s(£i-l)- ( 67 ) 

li=l,2 m=l s=0 i=2 

This establishes the one- dimensional character of i^HETC- The self-dualities 
of that Hamiltonian follow from those of the Ising ladder. 

It will be useful in the following to have an explicit description of the 
duality mapping connecting the two models of Eqs. (65, 67). First, we 



-"hetc 
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number the sites of the honeycomb lattice as in Fig. 3. The duality mapping 
then reads 



Wi 



(M,m) 



z z 
a ii,<f>(m) (T fi,<p(m)+li 



(68) 



so that $d(-^HETc) = -^h etc > see Fig- 2. The theory of operator algebras 
guarantees that isomorphisms like $ d of Eq. (68) must be equivalent to a 
unitary transformation [1]. It follows that there is a unitary transformation 
U d such that H° ETC = $ d (#HETc) = U d H HETC U ] d . 
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Figure 3: Numbering function (j> used to describe the duality between Hhetc and a one- 
dimensional model, illustrated for a lattice with L\ = 4 and L2 > 4. (f> numbers the 
hexagons, or equivalently, the sites of the underlying triangular Bravais lattice (large 
numbers). Since there is one point of the basis per hexagon, <j) numbers these as well 
(small numbers). The green crosses represent the transverse- field bonds a z present in 
^hetc- Periodic boundary conditions are considered. 

The Hamiltonian -£/hetc displays d = 1-GLSs that map into symmetries 
of the dual ladder model. Consider for concreteness the model illustrated in 
Fig. 2. It is clear that the operator 



(69) 



(using the numbering of spins described in Fig. 2) commutes with -^hetc- 
To determine the dual symmetry operator ^d^afa^c^) we need to compute 
$ d (crf ), which can be done according to the techniques of Ref. [1]. The result 
reads 

$ d K)=< 2 a 2 VM, 3 , (70) 
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so that the dual symmetry reduces to 



(71) 



The calculation of the last paragraph allows us to understand the restric- 
tions on the transverse-field bonds a z allowed in H^etc- That is, the second 
term in Eq. (65) included a sum over only one triangular sublattice. The 
introduction of omitted bonds such as erf on the other sublattice would de- 
stroy the exact dimensional reduction. To see this, notice that for a model 
with Li,L 2 arbitrary, the duality maps cr£ iai+a2+ t, to 



(Equation (70) is a special case of this result, and similar expressions hold 



for $dK iai+sa2+6 ), s = 2, ■ ■ • , L 2 ). So, not only the presence of a z Liai+a2+b 



in H\-\etc would introduce long-range interactions in H^ ETC , but the range 
of these interactions would grow without bound as a function of the area on 
which H\-\etc operates. This would spoil the exact dimensional reduction in 
the thermodynamic limit. In contrast, adding these bonds does not affect 
the d= 1-GLSs of the model. This shows that d-GLSs can exist even in 
the absence of exact dimensional reduction, as already anticipated in the 
introduction. 

The same crucial feature can be seen from a different perspective. It is 
intuitively clear that for any D-dimensional model one can always devise ar- 
tificially a dual, lower- dimensional, model. What sets most of these dualities 
apart from true exact dimensional reduction is again that the range of the in- 
teractions in the lower-dimensional model will in general grow without bound 
as a function of the size of the higher-dimensional model. Only dualities that 
fail to show this feature establish exact dimensional reduction. 

Thus, we conclude this section by highlighting its important result: the 
mapping of the D = 2 quantum system of Eq. (65) on a honeycomb lattice, 
with a staggered field appearing only on one triangular sublattice, into the 
d = 1 = d (ladder) short-range system of Eq. (67). 

6.2. Topological color codes 

We next briefly remark on dimensional reduction in the case of color codes 
[49]. The discussion will be very similar to that in previous section. 




(72) 



i=2 
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Color codes are defined on planar graphs such that three and only three 
links meet at any site of the graph (trivalent graphs), and such that any 
plaquette is bounded by an even number of links. There are two stabilizers 
per plaquette p: = Yliep a t with a = x or z. Thus, the operators B*' z are 
products of spins S — 1/2 at lattice sites that belong to a given plaquette. 
The total Hamiltonian is given by a sum of these operators (bonds) over all 
plaquettes 

htcc = -J2 b p-J2 b p- ( 73 ) 

p p 

We will not repeat the calculations explicitly for this model, they are anal- 
ogous to those carried out for the model investigated in Section 6.1 and the 
D = 3 model to be investigated in the next section, Section 6.3. Since, by 
definition of the lattice supporting a color code, all plaquettes share an even 
number of common sites [49], the bonds B° and B^ t commute always, 

[B°,B$] = 0, (74) 

for all p,p',a, and a', and 

(#p) 2 = I- (75) 

If the lattice is realized on the surface of a torus, the bond algebra is sup- 
plemented by two additional relations (aside from Eqs. (74, 75)). These are 
given by 

U b p =1 ( 76 ) 

p 

for a = x, and z. That is, the product of all operators B* over all plaquettes 
p is unity (and the same holds true for the product of all stabilizer operators 
Bp). The system of Eqs. (73, 74, 75, 76) is defined on a Hilbert space of 
dimension 2^ with N s the number of sites of the system. The problem, as 
formulated, has a bond algebra isomorphic to that of two decoupled, periodic 
Ising chains in zero transverse field (see Ref. [2] for an analogous duality for 
the toric code model). Thus, this system has an identical spectrum to two 
independent, periodic classical Ising chains. 

Color codes have d= 1-GLSs. In particular, the toric color code just 
described commutes with any string operator obtained by multiplying all 
the spins a x or all the spins a z on the links of any closed curve [49] . 
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6.3. The XXYYZZ model 

Our third example of dimensional reduction entails a system that we call 
the XXYYZZ model [50]. It can be described in terms of S — 1/2 spins 
placed on sites m of an fee Bravais lattice. For concreteness, we can describe 
the sites m = m 1 a 1 + 777 2 a 2 + m 3 a 3 , G Z with /i = 1,2,3, as linear 
combinations of the basis vectors, as illustrated in Fig. 4, 

e 2 + e 3 e 3 + e. x e x + e 2 , N 

«1 = o > a 2 = o , «3 = o • (77) 




Figure 4: The basis vectors ai, a,2, as for the fee Bravais lattice. For any site m in an fee 
lattice, there is one and only one octahedron with m as its lowest vertex. 

The distinguishing feature of the fee lattice that we would like to exploit 
is the fact that it has exactly one elementary octahedron per lattice site. In 
other words, for any site m of the fee lattice, there is exactly one elementary 
octahedron having m as, say, its lowest vertex (and with every other vertex 
in the lattice as well). Then, we can introduce the operator 

f\ x x y y z z (7H.\ 

— a rn+a 1 -a 2 a m+a 3 a m a m+e2 a rn+a 3 -a 2 a m+a 1 i \'°) 

associated with each site in the lattice (notice that e!,e2,e 3 are integral 
combinations of a 1; a 2 ,a 3 , and so they represent displacements in the fee 
lattice, even though they do not form a basis). In other words, the interac- 
tion terms are given by the product of all sites surrounding a given site of the 
even sublattice with the product being of the form "XXYYZZ" wherein the 
component of the spin appearing in the product is determined by its location 
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relative to the center of the octahedron formed by the six sites. This con- 
stitutes a, D = 3 generalization of the "XXYY" product appearing in Wen's 
D = 2 model [35] (that is also dual to Kitaev's Toric code model [19]). 

Next, we introduce periodic boundary conditions in all three directions 
ai, 02,^3, to obtain a finite fee lattice Af cc with the topology of a three- 
dimensional torus. More specifically, we identify any two points m and m' 
equal up to a multiple of the period L M in the a p direction. The fee structure 
forces the periods L M to be even integers as should become clear below (see 
Fig. 5). The Hamiltonian of the XXYYZZ model reads 

#XXYYZZ = -J °m- (79) 

meA£ c 

This Hamiltonian acts on a state space of even dimension 2^, with N s = 
L1L2L3 an integer multiple of four, and displays d= 1-GLSs that are de- 
scribed in Ref. [50]. Our goal is to employ the general techniques introduced 
in [1] to look for an exact duality of -/^xxyyzz t° a lower- dimensional model. 
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Figure 5: The fee Bravais lattice can be built as an ABAB stack of two square lattices 
(A and B), displaced relative to each other (we show a lattice elongated in the stacking 
direction to make the figure more clear). Since the planar lattices A and B are both 
bipartite, the sites of an fee lattice can be grouped into four disjoint classes. Any two 
octahedra belonging to sites in the same class share at most one vertex. 

The XXYYZZ model shares several features of typical models displaying 
TQO: its bonds O m commute with each other and the degeneracy of the 
ground energy level depends on the topology of the underlying lattice (that 
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we have fixed to be toroidal here). The bonds {O m \ m £ A^ c } generate a 
commutative bond algebra. We saw already that dualities are encoded in the 
structure of this bond algebra, so we need to specify it completely. The extra 
information we need is that, by virtue of the periodic boundary conditions, 
the bonds O m satisfy four constraints. 

To understand this, notice that we can split the sites of an fee lattice into 
four disjoint classes (at least for periodic and open boundary conditions). 
The reason is that, as shown in Fig. 5, the fee lattice can be constructed 
as a stack ABAB of bipartite square lattices A, B displaced relative to each 
other. Thus we can distinguish four types of sites, according to whether they 
belong to the sublattice A\ or A 2 of an A plane, or B\ or B 2 of a B plane (see 
Fig. 5). On the other hand, the bonds O mK in any one class k = A±, A 2 , Bi, 
or B 2 share exactly one site with its neighbors in the same class (and the 
corresponding spin, see Eq. (78)). It follows that 

II °m K = 1- (80) 

These are the four constraints that complete the specifications of the com- 
mutative bond algebra for Hxxyyzz- 

As explained in Sec. 6.1, a Hamiltonian will be dual to ^xxyyzz if the 
two models share bond algebras with identical structures (isomorphic bond 
algebras), and act on state spaces of the same dimension. Let us consider 
then the model 

4 N s /4 

HuP = ~J^2 ^2 a K,m aZ K,m+li ( 81 ) 
k=1 m=l 

that represents four independent (decoupled), periodic Ising chains. The 
bond algebra generated by the set of bonds {a z K m a z K m+1 \ k = 1, 2, 3, 4; m = 
1, ■ ■ • , N s /4} is characterized by the fact that all bonds commute, and satisfy 
four independent constraints 

Ns/4 

II <m< m+1 = 1, V«. (82) 

m=l 

Also, H^p acts on a state space of dimension (2 7Vs / 4 ) 4 = 2 Ns . It follows 
that i?4ip is dual (unitarily equivalent) to Hxxyyzz- To establish a duality 
mapping, let k = A±, A 2 , B ly B 2 correspond to k = 1,2,3,4, and let K be 
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any numbering of the sites m K in the sublattice k = Ai, A 2 , Bi, B 2 , from one 
to N s /A. Then 

Om K > (T K,</. K (m K ) <T K,0 K (m K )+l (^) 

is a duality isomorphism that maps -£/xxyyzz to H 4 \p, i.e., $d(-f^xxYYZz) = 
Wd-ffxxYYZz^j = i/4ip. It follows that the XXYYZZ model must display exact 
one- dimensional behavior. Thus, in conclusion, we mapped the D = 3 system 
of Eqs. (78, 79) on an fee lattice to the d = 1 = d system of Eq. (81) 
representing four decoupled Ising chains. 

It is important to stress that the source of the dimensional reduction in 
these two last models is not the completely commutative structure of their 
bond algebras. For instance, a model as simple as the D = 2 Ising model in 
zero magnetic field 

fM=l,2 r 

with r = r 1 e 1 + r 2 e 2 , e Z, denoting the sites of a simple square lat- 
tice, displays a completely commutative bond algebra but no dimensional 
reduction, due to the presence of constraints 

1 = o~ r o~ r+ei X o~ r+ei a r+ei+e2 X o~ r+ei+e2 or r+e2 X 0~ r+e2 0~ r . 

Notice also that the Ising model has no d-GLSs with d < 2. 

6. 4- Effective boundary theories 

In the exact dimensional reduction mappings that we studied thus far in 
this section, the dimension of the Hilbert space was preserved. Thus, for 
instance, when mapping a two-dimensional system onto a one- dimensional 
one, the number of lattice sites was unaltered. By contrast, in Section 3, 
we illustrated how local boundary theories can be constructed wherein the 
system size was reduced; these local boundary theories, however, only pro- 
vided bounds on the correlation functions and, generally, did not constitute 
exact duals to the original high-dimensional theory. In this brief subsection, 
we wish to briefly point out that in some cases local boundary theories can 
constitute exact duals to the bulk system. In such instances, the general 
approach of exact dualities that we considered earlier in Section 6 coincides 
with that of the boundary theories of Section 3. 
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Towards that end, we seek an effective Hamiltonian H e ff such that 



Z = Tr A (exp(-/3#)) = Af Tr r (exp(-/3iJ efF )). 



(85) 



Here, H denotes the full Hamiltonian of the D- dimensional system A = TUA, 
whereas if e ff denotes an effective Hamiltonian on a <i-dimensional boundary 
T. An overall normalization factor is indicated by M . 

In the quantum spin lattice models that we considered above, the trace on 
any set of chosen bulk fields (or bonds appearing in the Hamiltonian [18, 1]) 
can be done trivially. In the case of the system of Eq. (65) in the absence 
of an external field (h = 0) as well as the systems of Sections 6.3, and 6.2, a 
partial trace over some of the bonds renders the remaining system to be that 
of decoupled bonds (or Ising spins in this case). In particular, the effective 
boundary theories that contain all bonds that have their support on T alone 
is given by that of decoupled Ising spins, 



That is, the trace over the bulk fields A gives rise to a Hilbert space volume 
factor alone that multiplies the partition function associated with H e ff of Eq. 
(86). Now, the theory of Eq. (86) were to result from the system of Eq. (65) 
in the absence of an external field (h = 0) and the systems of Sections 6.3 
and 6.2, if we were to "freeze" out the bulk spins in a spirit similar to that 
the derived local bounding boundary theories of Section 3. 

7. Consequences of dimensional reduction: limits on the storage of 
quantum information 

In general, we may apply the bounds of [2, 3] to arrive at both upper and 
lower bounds on all correlation lengths and autocorrelation times in general 
systems with d-GLSs. In particular, these bounds on autocorrelation times 
provide limitations on quantum memory schemes and illustrate that, in some 
cases such as that of Kitaev's toric code system [19] , information will be lost 
after a period that does not scale with the system size. 

Kitaev's toric code model is a D = 2 system defined on a torus (as the 
name suggests) that exhibits a d = 1 Z2 x Z2 symmetry. The degrees of 
freedom in this model are spin 5 = 1/2 operators (afj or afj) on the links 
(ij) of a square lattice having | A| = N s vertices that is endowed with periodic 




(86) 
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boundary conditions. (To conform with the standard conventions, we denote 
the links by (ij) rather than the particular dual lattice site on which they 
are centered.) Specifically, this model [19] is defined by the Hamiltonian 

N s N s 

#Kitaev = — A s ~ B p , (87) 

s=l p=l 

with the star and plaquette operators A s and B p given by 

x n^j n 4- ( 8§ ) 

J (*i)ep 

In the first of Eqs. (88), j denotes all nearest neighbors of s. In Kitaev's 
model, information is stored in two "topological qubits" determined by the 
eigenvectors of the d = 1 symmetry operators 

z h2 = n 4, x li2 = n 0%, (89) 

where C^C^) are horizontal and C2(C[) vertical closed (d = 1) contours on 
the lattice (dual lattice) on the torus. As illustrated in [2], Kitaev's model 
is special in that it may be exactly mapped to two decoupled classical Ising 
chains (i.e., a system of dimension d = d = 1). Such a mapping is performed 
along lines similar to those of Section 6. The operators A s and B p map onto 
bonds of the two decoupled Ising chains (each of length N s ). The mapping 
is achieved by replacing, in the sum of Eq. (87), A s — > a z s a z s+1 on one chain 
and B p — > t z t z +1 on the other chain. That is, the Hamiltonian of the D = 2 
dimensional system of Eq. (87) is dual to that of two Ising chains (d = d = 1) 

N s N 3 

Rising-chain = ~ ^ ~ T P T P+l> ( 90 ) 

s=l p=l 

with o~ z N+l = o~\ and r^ +1 = rf . Similarly, within this mapping, the topolog- 
ical qubits Zi and Z 2 map onto single spins (say those located at the sites 
number 1 on both chains), Z\ — > a\ and Z 2 — > r*. The above set of mappings 
preserves all of the algebraic relations amongst the bonds {A s } and {B p } that 
appear in the Hamiltonian of Eq. (87) as well as the two qubits Z\ and Z 2 - 
We now discuss consequences of dimensional reduction on the storage of 
quantum information. First, we very briefly discuss how, physically, infor- 
mation may be lost by defects in this particular example system of Kitaev's 
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model. Then, we turn to a formal framework that illustrates how bounds on 
the storage of quantum information follow from our EQDR theorem. Lastly, 
we remark how such dimensional reductions allow for an understanding of 
particular cross-over temperatures in finite size systems. 

Conceptually understanding how dimensional reduction, when it occurs, 
may limit the storage of quantum information is relatively straightforward 
when thought about in terms of defects that eradicate the character of an 
initial stored state and lead to a finite correlation length. In Kitaev's toric 
code model, dimensional reduction implies that topological defects that arise 
in one-dimensional Ising systems rear their head anew on the D = 2 torus. 
Similar considerations apply for the other systems discussed in Section 6. 
In an Ising chain the proliferation of domain walls leads, at any positive 
temperature, to a finite correlation length and to finite autocorrelation (or 
memory) times. In accord with the one-dimensional character, the low tem- 
perature autocorrelation times may scale as r ~ exp[/3A] where A is an 
energy penalty associated with a domain wall. In a one- dimensional Ising 
chain with exchange constant J, this penalty is system size (N s ) indepen- 
dent and set solely by the energy scale J. (In Eq. (90), the nearest-neighbor 
exchange constant J = 1.) Thus, any information that is initially stored in 
Kitaev's toric code model will, at sufficiently long times, be destroyed by 
thermal fluctuations [2]. 

We now step back to consider the more general consequences of effective 
dimensional reductions borne by symmetries (Section 3). (In Kitaev's toric 
code model that exhibits the d = 1 GLSs of Eqs. (89) as well as other 
systems in Section 6, these effective dimensional reductions become exact.) 
To obtain bounds on memory times, we set the operator function / of Section 
3 to represent an autocorrelation function. That is, we consider what occurs 
when we choose / = Z 1 (0)Z 1 (t) so that its expectation value will measure 
the autocorrelation of an initial quantity Z\ with itself at later times t. If 
quantum memories are stable then (f) D does not decay at long times. By 
the results of Section 3, the autocorrelation function in a general system with 
d-GLSs, is bounded (both from above and from below) by autocorrelation 
functions of a <i-dimensional system in which the range of the interactions 
and symmetries are preserved 

<$L < G{2) (t) = <Si(0)£i(f)> < G( 2 p ) per , (91) 

(2) (2) 

with G\ a ^ er , Crupper autocorrelation functions in a d- dimensional system. These 
lower-dimensional autocorrelation functions serve as upper and lower bounds 
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on the autocorrelation function in D spatial dimensions. It is notable that 
in systems with interactions of finite range and strength the autocorrelation 
times of quantities not invariant under d = 1 GLSs do not scale with the 
system size. Thus, in such systems, the autocorrelation times are system size 
independent. In a similar fashion, system size dependent autocorrelation 
times appear in the ordered phase of low- dimensional systems with long- 
range interactions exhibiting a phase transition, such as ferromagnetic spin 
chains with l/r A interactions (A < 2) [53]. Taken together with the bounds 
of Section 3, this allows for divergent autocorrelation times of quantities that 
are not invariant under d-GLS in systems having long-range interactions. It 
is notable that even though excitations may involve a macroscopic number of 
spatially local degrees of freedom in systems such as Kitaev and the XXYYZZ 
model, and thus it might be expected that the system exhibits glass-like 
characteristics [50], what matters in the computation of the autocorrelation 
functions is not how memory may be erased in real space. Rather, what is of 
importance, so long as the noise is random, is the algebra underlying these 
excitations. Thus, we would like to suggest that, at least for some realizations 
of noise stemming from an external heat bath, the XXYYZZ model exhibits 
one- dimensional memory type effects and not glass-like features. 

A peculiar feature of Kitaev's toric code model is the existence of a cross- 
over temperature which scales with the system size [2, 54]. In what follows, 
we discuss this observation through the prism of the exact dimensional re- 
duction to Ising chains that Kitaev's model displays. In an Ising chain with 
nearest-neighbor interactions (of strength J = 1) at an inverse temperature 
0, the correlation length scales as £ ~ — 1/ ln(tanh/3). This scaling implies 
the existence of a crossover temperature scale T cross below which the entire 
chain is correlated and above which the chain length is larger than the cor- 
relation length. On a single Ising chain with N s ^> 1 spins, this gives rise to 
a crossover temperature T cross ~ 1 / In N s [2] . A similar crossover tempera- 
ture is seen when the entanglement entropy is analyzed [54] . A rudimentary 
understanding of this scaling can be seen from standard energy-entropy bal- 
ance considerations in the free energy, F = E — TS, for the insertion of a 
d = 1 defect (i.e., a domain wall in a one-dimensional chain for discrete 
d = 1 symmetries) in a system with short-range interactions. Associated 
with the insertion of a defect there is an energy penalty E = A and an en- 
tropy S ~ lniV,,. These two contributions lead to a a crossover temperature 
scaling as T cross ~ l/lniV s . For T < T cross defects are unfavorable, while 
T > T cross defects proliferate. The discussion above can replicated mutatis 
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mutandis for the systems of Section 6 when these reduce to Ising chains as 
well as many others which reduce to other low-dimensional systems. The 
existence of predicted crossover temperatures in these systems, as adduced 
from our general considerations above regarding dimensional reductions, may 
hopefully be bolstered by explicit computations of thermodynamic quantities 
[2] or entanglement entropy [54] in future works. 

8. Conclusions 

In conclusion, we examined both (i) effective and (ii) exact dimensional 
reductions and holographic correspondences. Specifically, we proved that 

(i) In a rather general local D-dimensional quantum theory that is defined 
on some volume A, correlation functions involving fields that lie on a subspace 
r C A of dimension d < D can generally be bounded by correlation functions 
in a local d- dimensional theory. These bounds lead to an effective dimensional 
reduction which become most potent when the original D-dimensional theory 
exhibits a particular symmetry (or set of symmetries) having its support on 
the region T. (A symmetry operator that has its support on a general region 
of dimension d is termed a d-dimensional gauge- like symmetry.) When the 
symmetry operator acts on the d = d spatial region T, the resulting bounds 
of Eq. (35) lead to a generalized quantum Elitzur theorem. 

(ii) Stronger than effective dimensional reductions are exact dimensional 
reductions that relate systems in different dimensions. We studied and illus- 
trated exact dimensional reductions of two different types: 

(a) Dualities centered on the large-n limit of 0(n) vector theories and the 
high temperature (or weak coupling) limit. The dualities derived in this case 
relied on the dependence of the partition function (or generating functional) 
on a local effective density (in the latter case, the density of states for sin- 
gle modes in the non-interacting limit which retains its character also in the 
presence of interactions in the large-n limit or that of high temperatures). 
We briefly speculated that it may be possible to search for analogues of func- 
tional theory type treatments [36] in field theories and interacting many-body 
systems that effectively render the theory nearly non-interacting. With the 
aid of these functionals, it may be possible to more generally construct trans- 
formations that do not deform the appropriate functions and thus preserve 
the features of a local theory when it is related to another local system in a 
different spatial dimension. 

(b) We further examined exact dualities that do not appear in some limit but 
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that, rather, always hold true. The studied dimensional reduction dualities 
focused solely on the physical bonds that appear in the Hamiltonian or action 
defining the theories. With the aid of this approach, redundant degrees of 
freedom on both sides of a general duality (e.g., gauge symmetries in various 
gauge theories, local coordinate invariance in gravity theories, etc.) may be 
automatically discarded. The equivalences in case (ii) focus on the form of 
the spectrum of the theory: case (a) focuses on the local density of states (in, 
for instance, fc-space) and (b) on mappings that trivially preserve the exact 
spectrum of the theory. 

Both the (i) effective and (ii) exact dimensional reductions introduced in 
this work do not require an actual physical compactification of dimensions. 
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Appendix A. Entanglement-based version of the EQDR theorem 

We briefly describe a variant of the EQDR theorem that holds for arbi- 
trary states p, at least if the Hilbert space of the system is finite-dimensional. 
The assumptions are again that we can treat the bulk and the boundary as 
distinguishable, Ha = %r ® Hx, and that operators / localized on Y have 
the form / = f-p £g> 1 A . Then we can expand an arbitrary state operator p as 



where pn, Pm are state operators, and Aj may be any real number, positive 
or negative such that = 1. Separable or unentangled states corre- 

spond to the case where all Aj happen to be positive. An entangled state is 
characterized by the presence of at least one negative coefficient in (A.l). 




(A.l) 
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Let us first justify Eq. (A.l). Notice that an arbitrary Hermitian operator 
can be written as the difference of two positive operators, by virtue of the 
spectral decomposition theorem. If O is Hermitian and Oj, Pi are its real 
eigenvalues and eigenprojectors, then 

O = °i p i = Yl 0i + Pi + ~ £ \°*- \ Pi - = 0+ - 0~, (A.2) 

i i+ i- 

where the index i + varies over that subset of i's for which Oj is positive 
(negative). We can identify the real space of Hermitian linear transformations 
Ja on Ha with the space Jr®J\, where Jr and J A , represent transformations 
on the boundary and bulk, respectively. Since these are finite dimensional 
spaces, it follows [31] that any Hermitian operator on "Ha, and any state 
operator p in particular, can be written as a real linear combination 

p = J20 ri ®0 Al . (A.3) 

i 

While there is no guarantee that the On, Ai will be state operators, we can 
apply the decomposition of Eq. (A.2) to each On, Ai to turn the expansion 
of Eq. (A.3) into an expansion in terms of products of positive operators. 
If we next normalize these positive operators, we get the expansion of Eq. 
(A.l), and we see that the Aj, although real, need not be positive by virtue 
of the minus sign in Eq. (A.2). 

We can now use Eq. (A.l) to compute inequalities of interest to dimen- 
sional reduction and holographies. If / = /r <8> 1 A is an observable 

L + (f)t - L_(/)f < Tr A (p/) < L + {f)t - L_</)-. (A-4) 
where L + = A i+ , L_ = |Aj_| are both positive, 

(f)u = max i+ Tr r (pri + /r), (/)" = min^Trp^p^/r), (A.5) 

and 

(/)+ = min i+ Tr r (pri + /r), (f)T = max 4 _Trp(p ri _/r). (A.6) 

Notice that if the state p is unentangled, then L_ = and L + = 1, and 
this version of the EQDR theorem, Eq. (A. 4), comes remarkably close to 
the classical dimensional reduction theorem, as perhaps is to be expected. 
In general, the inequalities of Eq. (A. 4) clearly illustrate another way of 
achieving effective dimensional reduction. 
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the classical covariant Hamiltonian as defined from the action vanishes 
identically. By contrast, the non-covariant Hamiltonian defined via a 
split of coordinates on a time slice can at times lead to a sensible results 
as in some early works on canonical quantum gravity (e.g., the ADM 
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heat, rendering the standard rules about stability of fluctuations sus- 
pect. For example, when energy is removed from a star, it gets hotter. 
The more massive a black hole is, the cooler it is, hence the Hawk- 
ing radiation is totally irrelevant for astrophysical black holes. It is for 
this reason that a black hole can never be in thermal equilibrium with 
an infinite heat bath, only with one of finite size. Thus, a number of 
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